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The Balitsky-Kovchegov (BK) evolution equation in its resummed in- 
tegral form as obtained in [U [2] is considered. We solve it numerically and 
compare to the unresummed BK equation formulated as an integral equa- 
tion and to the solution obtained by the BKsolver package. Sensitivity of 
the solution to an introduced resolution parameter and initial conditions is 
investigated. 
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1. Introduction 



Quantum Chromodynamics (QCD) is a theory which is used to set up 
the initial conditions for the collisions at the Large Hadron Collider (LHC) 
as well as to calculate properties of hadronic observables. At high energies as 
available at the LHC one enters into a region of the phase space where both 
the energy and momentum transfers are high and partons eventually form 
a dense system which is expected to saturate [3l IH El [7] • Indeed, there 
is gowing evidence that the saturation really takes place [H El El E]- The 
basic perturbative QCD equation which sums up the terms proportional to 
a n hi m (s/so) and also accounts for formation of the dense sytem of partons 
is the Balitsky-Kovchegov (BK) equation 0[9]. The BK equation, valid in 
asymptotic regime, does not take into account coherence effects in emission 
of gluons. This property manifests itself as independence on the hard scale 
associated with the external hard probe. Recently a framework has been 
provided in [H [21 [TO] [T2| [13] where both the dense systems and the hard 
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processes at high energies can be studied. This framework was based on the 
observation that the BK equation can be rewritten in an exclusive form and 
further extended to include the coherence effects. In the study presented 
here, which is a step towards understanding properties of the equations ob- 
tained in [H [21 [TO] , we perform a numerical analysis of a new form of the 
BK equation and compare it to the original formulation (see also [llj). In 
particular, we study the dependence of the resummed BK equation on a 
new scale which has the meaning of a resolution scale. We show that the re- 
summed BK equation agrees with the original one when the resolution scale 
li is already of the order of 10~ 3 GeV. We also investigate the dependence of 
the solution on the form of the initial conditions. 
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Fig. 1. The diagram explaining the meaning of the variables in the BK equation. 



2. Exclusive form of the BK equation 



At the leading order in ln(l/x) the BK equation [HJ [9] for the gluon 
number density in the momentum space is written as an integral equation 
and reads 00 [TO]: 
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where the lengths of transverse vectors lying in the transversal plane to the 
collision axis are k = |k|, I = |1| (k is a vector sum of transversal momenta of 
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emitted gluons during evolution), z = x/x'(see Fig. [TJ), a s = N c oi s /n. The 
strength of the nonlinear term is controlled by the targets radius R. The 
linear term in eq. ([1]) can be linked to the process of creation of gluons while 
the nonlinear term can be linked to fusion of gluons. The interplay of these 
two terms as a net effect leads to saturation of gluons. The unintegrated 
gluon density obeying the high energy factorization theorem |14) is obtained 
from [HI fl6]: 



J"bk(x, k 2 



■k 2 V 2 k $(x,k 2 ), 



(2) 



where the angle-independent Laplace operator is given by V? 
As shown in [H [2] this equation can be rewritten in a resummed form: 
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1 — k and Ar(z, k, //) = exp (— a s In i In is the Regge form 



Eq. is a form of the BK equation in which the resummed terms 
resulting in the Regge form factor are the same for the linear and nonlin- 
ear part. This form served as a guiding equation to generalize the CCFM 
equation |17| \TE[ I19| to the KGB JS equation [TJ [2] which includes nonlinear 
effects. These effects allow for recombination of partons with a constraint 
on an emission angle. 



3. Computational method 

In order to study the behaviour of both the equations Q and Q, they 
were solved numerically following an iterative procedure, which we detail in 
this section. 

For numerical treatment of Eq. ([Tj), the inner integral is approximated by 
reducing it to a finite interval (ko,qf). This allows for a direct numeri- 
cal representation of the function, which is the solution to be computed. 
Iterative refinement of it is more naturally expressed with w = — as the inte- 
gration variable. Also, to give an accurate view of the numerical procedure, 
we explicitely denote the second argument of $ as &2 = k 2 , 1% = I 2 , so the 
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equation reads 

l 2 ®(w,l 2 ) - k 2 $(w,k 2 ) k 2 $(w,k 2 ) 

\k 2 -i 2 \ + y/±q + k\ 

a s f xo dw „ , , , lS 
7rR z J x w 

where for simplicity we assume R = ^= and a s = Nc " s with the fixed QCD 

coupling constant a s = 0.2. The evolution starts at xq = 1CP 2 . 
The function $ is represented on a regular mesh of points (x m , k 2n ) that are 
equidistant in a logarithmic setting: logx m = log x m i n +mA x (and simlilarly 
in ^2)- For values of <3? with arbitrary arguments we employ a bilinear inter- 
polation with respect to logarithms of the variables. The iterative procedure 
starts with $0 as a first approximation of the solution. Then the right hand 
side of the equation is calculated repeatedly to provide and use subsequent 
approximations. Since Eq. ^Q) can be cast as a differential equation in x as 
well, the above algorithm can be modified to be more efficient. Namely, the 
particular structure of the equation permits us to work with the domain of 
$ reduced in x to cover the arbitrarily small interval (x',Xo). In order to 
extend some known solution to reach some smaller x', it is enough that the 
iterative recalculation of the right-hand side is performed only for the newly 
considered values of x. Starting with x' = xq and going down to next grid 
points we save computational time. In our implementation, the integrals 
are computed using the VEGAS method, as available in the CUBA library 
|20| . To ease the integration, extra variable changes: dw/w — > d (Inw) and 
dl 2 /l 2 d (ln/2) are employed. 



+ a s 



X0 ^y, 

w 



q2 f dh 



The following two forms of the driving term $0 are considered here: 

^(x,k 2 )=e W (-k 2 /GeV 2 ), (5) 
$ b (x,k 2 ) = {k 2 /GeV 2 y 1/2 . (6) 



The resummed form of the BK equation, eq. ([3]), is solved in approxi- 
mated form for similar reasons. While the lower limit of the inner integral 
is now exactly fj,, the upper limit and the vector sum |k + q| = S(q 2 , k 2 ,tp) 
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need to be constrained: 

<f>(x,k 2 ) = ®o(x,k 2 ) 

f x ° dw f^f dq 2 f 71 

+ a s — A R (z,k 2 ) — -/ dip$ (w,m(S(q2,k 2 ,<p),ko,q 2 f )) 
Jx w Ju? Kq 2 J J 



a. f x ° dw 



' 1 —A R (z, k 2 )^ 2 (w, k 2 ). (7) 



ttR 2 

The length of the vector sum 



S(a 2 , b 2 , a) = a 2 + b 2 + 2y a 2 b 2 cos a (8) 



may fall outside the ytQ,q 2 J range, hence the integrand is approximated by 
limiting the second argument of This is done as follows: 



m (k' 2 ,ko,q 2 f ) 



/cq when k' 2 < k 2 , 



k' 2 when kl < k' 2 < q 2 , (9) 
q 2 when k' 2 > q 2 . 



This way the domains of $ and of the integrals all can be kept finite. 

The initial conditions resulting from the resummation procedure include 
an extra factor as follows: 

^ b (x,k 2 ) = X^^\x,k2), (10) 

where Lu(k 2 ) = a s In 

Except for the above remarks, both equations are solved using the same 
method. The numerical solutions presented below are computed with k$ = 
0.001 GeV and qf = 100 GeV, for the cut-off values ranging from ji = ko to 
y. = 0.01 GeV. 



4. Numerical solutions 

In this section we present results of numerical solution of the considered 
equations. In Fig. [2] we present the solutions of the BK equation formu- 
lated as an integral equation compared to the more commonly used integro- 
differential formulation. The former is solved via the iteration method, while 
the later by using the BKsolver package |21j . As can be seen, the solutions 
stay within 1% for k < 1 GeV. Above this value, the solutions diverge and 
at the lowest x considered the relative difference reaches 4%. For the largest 
values of k, the distribution ^ approaches zero anyway, thus the agreement 
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Fig. 2. Upper plots: the solutions of the BK equation as formulated in eq. |T|) 
compared with the solutions obtained with the BKsolver package |21| . Lower plots: 
ratios of these solutions. 



between the two solutions is satisfactory. As the next result we present in 
Fig. [3] the solution of the resummed form of the BK equation and compare 
it to the unresummed one. In the former, the scale /x introduces numeri- 
cally some weak dependence of the solution on its value. This dependence 
is expected formally to disappear in the limit fi — > 0. In Figs. U and [5] we 
study the dependence of the solution of the BFKL and BK equations on 
the parameter /i. As already mentioned above, the resummation procedure 
assumes that the scale \x is the smallest scale in the problem and that it 
should tend to zero. One can see that the BFKL equation is more sensitive 
to the resummation parameter than the BK one. The reason for this is the 
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Fig. 3. The comparison of the solutions of the BK equation as formulated in eq. |T|) 
with the solutions of its resummed form given in eq. ([3]); fi = ko = 0.001 GeV. 
Lower plots: ratios of these solutions with values obtained for eq. §3§ in the de- 
nominator. 



feature of the saturation scale which provides a cut-off on small momenta 
and therefore weakens the dependence on [i. Particularly interesting is the 
shape of the gluon density from the BFKL equation shown in the upper 
right-hand side panel of Fig. [JJ One can see that when the cut-off is larger 
than the probed kt value, the distribution bends upwards. This is due to 
the fact that the Regge form factor in that case becomes larger than one. 
When we lower the cut-off, this structure disappears. Finally, in Fig. [6] we 
investigate the dependence of the solutions on the form of the initial con- 
ditions. The comparison of the BFKL and BK cases indicates that the BK 
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Fig. 4. Sensitivity to the cut-off /i without the nonlinear term. Upper plots: the 
solutions of the resummed BFKL equation, i.e. eq. ([3]) with the nonlinear term 
neglected; their ratios are plotted on the lower panes. 

equation is more universal, i.e. the spread of its solutions is smaller. 



5. Summary 

In this paper we have performed the study of the resummed form of 
the BK and the BFKL evolution equations. We have compared the ob- 
tained solutions to the unresummed ones as well as to the solutions provided 
by the BKsolver package. The solutions of the BK equation formulated as 
the integro-differential equation and of the BK equation formulated as the 
double-integral equation agree well, with the discrepancy staying below few 
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Fig. 5. The effect of the nonlinear term on the sensitivity to the cut-off /i. The 
solutions of the resummed BK equation of eq. ([3]); for \i = 0.01 GeV (red and 
next lines) and /i = 0.001 GeV (blue and next lines). Lower plots: ratios of these 
solutions. 

percent. We notice the residual dependence of the BK equation on the re- 
summation parameter and also observe that it is less sensitive to the resum- 
mation parameter [i as compared to the BFKL equation. We attribute this 
phenomenon to the emergence of the saturation scale in the BK equation. 
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Fig. 6. The comparison of the solutions obtained with different forms of the 
$o term. Left: the resummed BFKL equation, i.e. eq. ([T]) with the non- 
linear term neglected; right: the resummed BK equation of eq. ©; red lines: 
$ = x°exp (-fc 2 /GeV 2 ), blue lines: $ = x° (k 2 /GeV 2 )- 1 / 2 . 
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